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Abstract 

The problem of derivation of the equations of motion from the field equations is considered. 
Einstein's field equations have a specific analytical form: They are linear in the second order 
derivatives and quadratic in the first order derivatives of the field variables. We utilize this 
particular form and propose a novel algorithm for the derivation of the equations of motion 
from the field equations. It is based on the condition of the balance between the singular terms 
of the field equation. We apply the algorithm to a nonlinear Lorentz invariant scalar field model. 
We show that it results in the Newton law of attraction between the singularities of the field 
moved on approximately geodesic curves. The algorithm is applicable to the A'^-body problem 
of the Lorentz invariant field equations. 

PACS: 04.25.-g, 45.50.Dd, 95.10.Ce 



1 Introduction 

General Relativity (GR) is unique among the class of field theories in the treatment of the equa- 
tions of motion. When the particles are modeled by the singularities of the field, their motion is 
completely determined by the field equation. By comparison, in classical electrodynamics, the equa- 
tions of motion of matter sources are postulated independently from the field equations. A precise 
analysis of the basis assumptions of classical electrodynamics was given recently in the premetric 
axiomatic framework of Hehl and Obukhov [2]. This way, two independent assumptions have 
to be postulated separately from the Maxwell field equations: (i) The agent of interaction — the 
Lorentz force expression, and (ii) The agent of inertia — the mass-times-acceleration term. 

In GR, the field equations are nonlinear so a superposition of solutions does not satisfy the 
field equations. Consequently, the motion of sources cannot be completely independent of the field 
equations. Moreover, it was shown already in the early days of GR that the motion of a point-like 
particles is determined by the field equations and should not be postulated separately. This result 
was achieved in two essentially different ways: 

(i) Due to Einstein and his collaborators |S| — 01 (see also ^Hl — the equations of motion 
of the point-like particles were derived from the vacuum field equations 

R^J~\g^JR=Q■ (1) 

Point-like particles are represented by the singularities of the gravitational field, i.e, of the solutions 
of d. 
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(ii) Fock ^21) see also ^7], considered particles as represented by a suitable energy-momentum 
tensor Ty . Correspondingly, the equation of motion was derived from the full Einstein field equation 

R^l - IdvR = T^j . (2) 

Papapetrou |18j , does not assume a specific form of the energy- momentum. Thus his approximation 
method is fairly general and, in principle, can be carried out to arbitrary high orders. 

In both models, a suitable approximation scheme was used so that the equations of motion of 
massive particles, were approximated by geodesic curves. 

The Einstein-Infeld-Hoffmann (EIH) equations were used for high order PPN approximations of 
the TV-body problem — for calculation of the gravitational radiation reaction force |23|. 
and for describing the motion of a gyroscope QU- For the geodesic line approximation of the 
motion of bodies of a sufficiently small size and mass, see |26j and the reference given therein. 

The EIH-approach was used to establish arguments for the plausibility of the geodesic postulate. 
In particular, Sternberg proved the geodesic postulate in a general form for a system that 
incorporates also Yang-Mills fields. 

The aim of our paper is different. In this article, we adopt Einstein's idea of describing point-like 
particles by the singularities of a field. Then, some natural questions can be considered: 

(i) What analytical properties of the Einstein field equations are used in order to derive the proper 
equations of motion for singularities? 

(ii) What parts of the field equation have to be identified with the agent of inertia and the agent 
of the force"? 

(iii) What other nonlinear field equations lend themselves to derivability of the equations of motion? 

In this article we submit partial answers. First we define a suitable class of the field equations. 
These are linear in the second order derivatives and quadratic in the first order derivatives of the 
field variables. We also require the equations to be Lorentz invariant. We show that the linear part 
serves as the agent of inertia while the quadratic part is the agent of the force. 

We suggest a nonlinear scalar Lorentz invariant field model. There is an exact static solution 
to this equation, which is singular at N points. This solution does not imply interaction between 
singularities, it is non physical. Indeed, for the Newtonian gravity as formulated by a field theory, 
A = 0, a similar exact solution^ (j) = '^^rrii/ri exists. It is rejected by an additional axiom, 
which is Newton's law of attraction. Also for Maxwell electrodynamics, the potentials satisfy (in 
the Lorentz gauge) the wave equation OA — 0. An exact static solution is A = Qi/i^i- Also this 
solution is non physical (it does not describe the field of free non-constrained sources) because it 
contradicts Coulomb law. In both examples, one cannot distinguish between the quasi-static and 
dynamic solutions by the field equations. 

Wc do not prc-postulate the law of interaction. It is derived from the field equation. For our 
model, the non-physical exact static solution is rejected, the same way the quasi-static solution was 
rejected in Newton's and Maxwell's theories. Instead, an approximate solution is constructed. When 
our algorithm is applied to this solution, it yields the Newton-type law of attraction. 

Our algorithm can be briefly described as follows. 

(i) Start with a nonlinear Lorentz invariant field equation. Let its static spherical symmetric 
solution (exact or approximate) be singular at a point. Due to Einstein's description, this is 
the field of one particle. 

(ii) The field of a particle moving with a constant velocity is constructed by a Lorentz transfor- 
mation of the static field. 

(iii) The field of a particle moving on a curved trajectory^ ''/'(O is constructed by a Lorentz trans- 
formation with an instantaneous velocity xp{t). 

^We use a system of units with c = 1 , G = 1. 
■^The bold symbols denote 3-dimensional vectors. 



(iv) The field of N particles is taken to be, approximately, a superposition of the fields of the single 
particles. 



(v) When this field of N particles is inserted into the field equation, two singular expressions 
emerge: one from the linear part and the other from the quadratic part. By equating the 
highest order terms of these expressions the strength of the singularity of the solution is 
reduced. The corresponding balance equation yields Newton's law of attraction between the 
singularities. 

We apply this algorithm to a nonlinear Lorcntz invariant scalar field model. In this context, the 
Newton law of attraction is derived. 

The paper is organized as followed: In section 2, we construct a nonlinear scalar Lorentz invariant 
field model. It shares some analytical properties of the Einstein field equation. In section 3, exact 
static solutions are presented. The solution with N static singularities does not allow to describe 
the interaction. We treat it as non-physical. In section 4, we construct a non-static solution with N 
singularities. In section 5, we show that the balance between the singular contributions yields, in the 
first order, the Newton law of attraction. In section 6, we present a qualitative technical summary 
of our algorithm. In section 7, we give a brief outlook of the proposed alternative algorithm and its 
possible extensions. 

2 A nonlinear scalar model 

2.1 Pointwise singularities in GR 

In vacuum, the Einstein field equation can be written as 



Here and in the sequel, the Greek indices refer to the coordinates, ^,v, - ■ ■ — 0, 1,2,3, while the 
Roman indices denote the specific singular point of the solutions, — I,-- - ,N. From an 

analytical point of view, is a system of ten equations for ten components of the metric tensor 
g^i,. All the equations are: 

(i) Hyperbolic nonlinear PDE of the second order; 

(ii) Linear in the second order derivatives of the metric tensor; 

(iii) Quadratic in its first order derivatives. 

The field of a pointwise source is described by the Schwarzschild solution of Q, which nonzero 
components are: 

(i) Smooth for r ^ 0; 

(ii) Singular at the origin r = 0; 

(iii) Have the Taylor expansion of the form 



The positive parameter m in Q is interpreted as a mass of a pointwise source located at the origin 
of the coordinates r = 0, while the constants Ci, C2, • • • are dimensionless. 

Due to the translation invariance of Q, the mass m can be located at an arbitrary point, r = Tq. 
Moreover, it was proved, see that the system has a smooth solution with N singularities 
located at the points r = ri{t). Although, even for = 2, an exact form of such a solution is not 
available, the leading terms of its Taylor expansion can be calculated. The main output of the EIH 
analysis, is the fact that the functions ri{t) describe the proper motion of small bodies. Particularly, 
in the first order approximation, the functions ri{t) satisfy the Newton equation of motion with 
the Newton force of attraction at the right hand side. The higher order approximations yield the 
geodesic equation. 



(3) 




(4) 



2.2 Nonlinear scalar model 



In order to clarify how the field equation itself can determine an interaction between singularities, we 
will study a scalar field model which shares the special analytical properties of the Einstein equation 
mentioned above. 

Let a scalar field tp = f{x) be given on a flat Minkowskian 4£'-space with a metric -q^i, ~ rj^" — 
diag{\, —1, —1, —1). We consider a nonlinear scalar field equation 

T]^" iv^.fj^i^ - kip^f^ip^^) ^0 , (5) 

which is equivalent to 

Dv-k{ip^-\\Vip\\^)^0. (6) 

Here fc is a dimensionless constant, it's value will be discussed in the sequel. 

The commas denote partial derivatives (p^^ = d(p/dx^, the wave operator is □ = d^/dP — 
/dx^ — /dy^ — d'^/dz'^, V is the spatial gradient operator, the scalar product < • • • > and the 
corresponding norm 1 1 • • • 1 1 are Euclidean. 

In our approach, the scalar model (jsj is used only as a simplified illustration to the way an 
interaction between singularities can emerge in GR. The model, however, is interesting by itself. In 
particular, the numerical evolution of a scalar field model similar to the one under consideration was 
recently studied in 



3 Exact static solution 
3.1 Linear equation 

In order to derive a suitable exact solution to 0, we start with a special case k = 0. For this value, 
the equation |SJl is linear 

□ ^ = 7?'^"^,^, = 0. (7) 
Its unique static, spherically symmetric, asymptotically zero solution is 

m , , 

^ = (8) 

Note that this function is smooth for r 7^ and singular at the origin of the coordinate system. It 
is natural to interpret such a solution as a field of a pointwise particle with a mass (or a charge) 
equal to m located at the origin. The translational invariance of Q), yields that the singularity can 
be located at an arbitrary point Tq 

m , , 

^ " i? ' 

where 

fi = r-ro, R^\\R\\. (10) 

The linearity of the field equation Q yields the existence of a static solution which is singular at N 
distinct points 

I— I 

where 

R,=T-r,, R,^\\R,\\. (12) 

This solution approaches zero at infinity and is approximately spherically symmetric in small neigh- 
borhoods of every singularity. For long distance from the set of the singular points, the full spherical 
symmetry is reinstated 

M ^ 
= — , M - ^ m, . (13) 

1=1 

We interpret Hll() as a field generated by a static configuration of N particles of masses which 
are located at the points r^. Since the solution 111|) is static, an interaction between the particles is 
absent. 



3.2 A nonlinear equation 



Let us derive now a corresponding static solutions of the nonlinear equation with arbitrary values 
of the parameter k. We redefine the scalar field as 

(^ = -^ln0, (l)^e-'"^. (14) 



k 



It follows that 



?r iv.^^. - k^.f^V..) = -— ^'^"0,^- = . (15) 



1 

Thus, under the transformation p4|l . the nonlinear equation (|SJ) is transformed to the linear one Q. 
Consequently, the field <p satisfies the equation Q if and only if the new field cj) satisfies Q. We use 
this transformational property to derive certain exact static solutions of Q. 
Start with a solution of Q of the form 

Tfl 

0=l-fc-, i?=||r-ro||. (16) 

Under H14|l . it transforms to a 1-singular solution of ^ 

.= -lln(l-«=^). (17) 

This solution vanishes at infinity as m/R. Also, in the limiting case fc ^ 0, the solution H17|) 
approaches the Newtonian potential. Although these limits hold independently of the sign of k, for 
fc < 0, (|17|) is singular only at a point r = Tq and smooth at other points. We will see in the sequel 
that precisely the negative values of the parameter k yield the attraction between the singularities. 

Due to the transformation (|14l) , the solution lll|l with N singular points also can be transformed 
to a static solution of We take it in the following form 

i—l 

where ili = r — while is a radius vector to the i-th point. Some analytical properties of this 
solution can be motivated from the physical point of view. It is singular at the points r = and 
smooth otherwise. In the small neighborhoods of the singular points, the solution is approximately 
spherically symmetric. This full spherical symmetry is reinstated on a long distance from the set of 
N singular points, where the solution is approximated by 

N 

= --ln(^l-/c— ) , M = ^m, . (19) 

1=1 

In other words, the masses of the sources are approximately additive. Since the solution (|18|l is 
static its singularities do not interact. 



4 Non-static solution generated 

4.1 Instantaneous dynamical Lorentz transformations 

The solutions exhibited above are static. In order to describe an interaction between singularities 
we are going to derive a certain set of dynamical solutions with the same type of singularities. For 
this we have to move the singular points on curved trajectories. The corresponding field if will also 
move. For a Lorentz invariant equation the dynamical solution can be derived by a corresponding 
Lorentz transformation. 



Let us start with a unique singular point which moves on a straight trajectory with a constant 
velocity v. This motion can be represented as a Lorentz transformation of a stationary point r with 
the velocity —v. Using the Lorentz transformation with the velocity — v we have (see Appendix A) 

r^r-Q;v<v,r> - j3wt , (20) 

where the Lorentz parameters are denoted by 

To represent a motion of N singular points on distinct straight trajectories, we apply, in their 
neighborhoods, N independent Lorentz transformations 

Ti - a^Vi < Vi,ri > -/3iV^t, (22) 

with the Lorentz parameters 

P^ = ' = ^ I 1 - I • (23) 

In order to construct a time dependent solution, we apply, for the static solutions, an instanta- 
neous Lorentz transformations with time dependent velocities. Every such transformation acts only 
in a neighborhood of the corresponding singular point of the mass rui . In particular, for a curved 
trajectory ■0(t), we substitute in (|22() v by i/) and wt by tp- Consequently, a dynamical Lorentz 
transformation which depends on the trajectory is 

T ^ V - aip < ip,v > -j3ip ■ (24) 

In order to deal with a motion of of N distinct singularities, we define a set of instantaneous 
dynamical Lorentz transformations (i = 1, • • • , N) 

(r - n) -> (r - r,) - aitp^ < tp^, (r - r,) > , (25) 

where the Lorentz parameters and Pi are functions only of the velocity 




4.2 Non-static solution of the linear field equation 

Now, consider the linear field equation 

^'^>M- = 0. (27) 

Recall that this equation is a special case, fc = 0, of our non-linear model ((S)). 
We start with an exact static solution of (|27|l 

N 

^ = E ^ ' R^^rr,. (28) 

i—l * 

The transformation H25|l leads to a time depended ansatz 

N 

'^ = 11^' = (r-r,)-a.V>, < V>„(r-r,) > -ft^,. (29) 



Clearly, for an arbitrary function H29|l does not satisfy H27I) . Let us look for functions ipi{t) so 

that (|29|l is a solution of 127() . at least approximately. Calculate, approximately, the left hand side 
of (|27|l for If given in (|29|l . In particular, we omit the time derivatives of the Lorentz functions ai 
and (3i, which bring only terms 0(xl)ip). In this lowest order approximation, cf. Appendix B, 

° ^ = E 1^ < > • (30) 

i—l * 

Consequently if is an approximate solution of the field equation if and only if i/'j ~ 0. Thus, 
the linear field equation {Tj) can serve as a field model for a free inertial motion of a system of N 
singularities. 

In fact, for ■0, = 0, i.e. for xp^ — xpjt, this result is exact. Indeed, at the point Ri = 

r - r,: = i/), < ip^, (r - r^) > +l3itj , (31) 

which results in the equation of motion of a free particle 

r-r,^iP,t. (32) 



4.3 Non-static solution of the non-linear field equation 

Let us now turn to the full non-linear scalar field equation 

V'"'{v^^..-k^,^.^,u)=0. (33) 
If the transformation H25I) is applied to its static solution 

1 ^ 

4l-(l-fcES' R.-r r.. (34) 



. -I Ri 



we will obtain no more than the inertial motion i/) = 0. This is because the equation l|33() is Lorentz 
invariant. Instead of that, we start with a corresponding approximate solution, 



N 



'P = -^El^M^"^^) ' R^^r-r,. (35) 



z— 1 ^ 

which is a superposition of N independent exact solutions. 

In order to generate a non-static (interaction) solution, we apply, in a neighborhood of every 
singular point, Lorentz transformations with variable velocities. Thus we come to 

<^ = -^Eli^(l-^9i) ' -R. = (r - r,) - a.V', < ^„ (r - r,) > -A'/'. • (36) 

As in the linear case, for an arbitrary function xp{t), this ansatz is not a solution of H33|) . Let us 
now search for conditions under which (|36|l turns into an approximate solution of (|33|) in the lowest 
order. 

Calculating for H3()|) the linear second order derivative part of (|33|l . we get (see Appendices C 
and D) 



< Ri , Rj > 



< i)„ Rj > < i/),, Ri > {(3il3j + a, + aj - a^a^ < ipi, -0 ■ >) + 0{ipi)) . (37) 



The expressions are calculated up to 0{xjjil)). To this accuracy, the derivatives of the Lorentz 
parameters ai,(3i are zero. 

The second line of H37|l is of the order 0{ipj^ipj). For the lowest approximation, we neglect this 
term. Moreover, to the same accuracy, we can take (3i — I. Thus, from this stage, our results will 
be valid only for slow motion of the singularities. 

To sum up, we derived that the ansatz ()36|l is an approximate solution of the field equation H33|) 
only if the functions tpiit) satisfy 

^m, < > _ , 1^ ^ n n -> dR-) 

l^El -"^l.T3I|IT3W<^-^^>' (38) 

i=i » Ri i^j Ri Rj 

where 

R, = (r - r,) - a,^, < V-,, (r - r,) > -ft'i/', . (39) 



5 An interaction is generated 

Both sides of the equation (|38|l are functions of a point r. It is singular at N points = and 
regular at other points. Let us examine how this relation can be satisfied in a neighborhood of a 
singular point. Consider the p-th singularity. We will examine the field at an arbitrary point r close 
to this singularity, i.e., we have Rp 0. For the i-th singularity 

Rp^O yields R^ ^ R^p for i p , (40) 

where Rip is a vector directed from the point i to the point p. 

On the left hand side of the equation H38|l . there is only one singular term 



TTlp < Rp,-4^p> 



For the limiting values, we have in H4(J|I 



(41) 



The singular term on the right hand side of H38|l is 

rdp rrij 

-^Tzf^ET^^<^p'^^ > ■ (42) 

Rp j^p Rj 

The terms and (|42|1 are of the order 0{R~^) near the j>-th singularity. 
They cancel each other if and only if 

rrij 

< Rp,i)p>^ -k^^-^^ < Rp,Rj > . (43) 

j^p Rj 

Take into account that the point r (in Rp) is arbitrary. Hence (|43l) is valid only if 

ITLj 



j^P Rj 



i>p = -ky ^ ^ Rjp . (45) 



For distances greater than the Schwarzschild radius r = km, (|45() remains in the form 

^p = -^E|r^.p- (46) 



For a system of two singular points, it yields 



, mi 7712 _ 

miV'i = fc-^ i2i2 . (47) 



For fc < 0, H46(l and H47(l result in attraction between the particles. The absolute value of k is 
unimportant, since it amounts to the rescaling of the mass. 

This way the Newton-type law of attraction is derived from the scalar field equation. 



6 Qualitative Description of the Algorithm 

In this section we give a qualitative description of the proposed algorithm for deriving the equations 
of motion from the field equation. Consider a field equation of a general type 

a$,^, - = (48) 

without specification of the tensorial nature of the field variable <I>. Here, the coefficients a and h 
are the dimensionless functions of the field variable (or constants) 

a = a($), 5 = 6($). 

We do require, however, that H48() is Lorentz invariant. 

Let the field equation H48() has a static spherical symmetric solution $(r — Tq) with a singularity 
located at r = Tq. Denote the Lorentz transformation based on the velocity v by i(v). Conse- 
quently, if $ = <f>(r — ro) is a time independent solution of H48|l then L(v)$ is also a solution of the 
same equation for an arbitrary Lorentz transformation L(w). Note that if $ is a multi-component 
quantity like a tensor then L(v) involves not only a coordinate change but, also, a transformation 
of components of The solution i(v)<i> describes the field of a pointwisc singularity moving with 
a constant velocity v on the trajectory ij} — vt. 

Let us try to construct a generalization of the Lorentz transformation used above such that the 
origin moves on a curved trajectory xp ^ xp{t). Denote such a transformation by N{xp). The choice 

is a plausible candidate. Correspondingly, N{xp)^ is a rigid motion of the field. 

Now substitute N{2p)^ in igHJi. U -ip = const then A^('0)$ is also a solution of If not, then 
the linear part produces extra terms. These extra terms come from two sources: 

(i) The derivatives of the Lorentz root y^l — Ixpl" . In our consideration the velocity of the particle 

ip as well as its time and spatial derivatives are assumed to be small. It follows that the 
derivatives of the root are of the form C'(|'0| |'0|). Thus they may be rejected. 

(ii) The linear part. Since is time independent the first order derivatives of $ are only the spatial 
ones. Thus, the first order derivatives of L{ijj)^{x) involve spatial derivatives of <I> multiplied 
by ip. The second order derivatives of L{ip)^{x) cancel each other exactly in the same way as 
in a Lorentz transformed solution L(v)$(a;). One exception: The second order derivative ^p 
multiplied by spatial first order derivatives of $ do remain. This extra term that comes from 
the linear part of the field equation is the agent of inertia. 

(iii) The quadratic part. It involves only first order derivatives of the field. Consequently the fact 
that is variable does not affect it's form. 

Construct now a solution that describes a field of N particles, i.e. a field with N singular points. 
It can be approximated by a superposition of 1-singular solutions moving on arbitrary trajectories 

<l> = ^L,(V,,)$,(a;). 
j 



Substituting this approximate solution in (|48|l we obtain, in the hnear part, only the second deriva- 
tives xpj multiplied by the spatial first derivatives of 

Consider the quadratic part. It is composed of the first order derivatives of Lj{ij}j)^j{x) multi- 
plied by the first order derivatives of Lk{'tpk)^k{x). If j = k, since ^k{x) is a solution of (|^ . these 
products are cancelled by the linear part operating on li j ^ k the products will be declared 

to be an approximation to the interaction between the j-th and the fc-th particles. 

Near the fc-th singularity, for the linear part, only the terms coming from $^(2;), similar to (|41|l . 
will be dominant. Likewise, for the quadratic part, only the terms involving the derivatives of $fc(a;), 
similar to (|42|l . will be dominant. Equating, near the singularity, of the two terms above (again to 
the leading order) should, hopefully, result in the Newton-type law of attraction. 

7 Conclusions 

The motivation for our paper is to learn what properties of the Einstein field equations lead to 
derivability the proper equations of motion for the singularities. We conjecture that it pertains 
to the specific structure of the field equations: As partial differential equations, they are linear 
in the second order derivatives and quadratic in the first order derivatives of the field variables. 
Furthermore, we conjecture that the quadratic terms are the agent of interaction, they have to 
generate the force expression. The second order derivatives play a role of the agent of inertia, it has 
to generate the mass-times-acceleration term. 

We curry out the above plan for a simple non-linear scalar model. It is remarkable that only with 
the slow motion assumption, one is ultimately led to the Newton-type law of attraction. It should 
be emphasized that our scheme does not pretend to be more precise than the alternative methods 
cited above. Instead, we are looking for a proper identification of the terms of the field equation 
with the inertial and interaction terms of the equations of motion. 

In order to apply the proposed scheme to Einstein's equation, one has to solve the problem of 
consistency . Every one of ten independent equations Rp,i, = has to give the same equation of 
motion. Also the Bianchy identities have to be taken into account and their role in our scheme has 
to be revealed. It is well known that the Bianchy identities play a crucial role in EIH-procedure. 

An other extension is in the direction of embedding |3(Jj . The trajectories obtained are approxi- 
mate. At this point, two avenues are open. The first one, which is adopted by EIH is to get higher 
order approximations to the trajectories. This procedure is also used in the PPN approach. By 
these methods, the successive approximations become highly singular near the particle trajectories. 
The second avenue is to embed the singularities in a field satisfying the field equations. For that 
purpose, the successive approximations should add regular terms (and, possible, low order singular 
terms) near the trajectories. 

The method proposed in this paper can hopefully be useful for a rigorous proof of the geodesic 
postulate. 
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A Lorentz transformations 

Lorentz transformations are usually written in a very special form when the axes of two reference 
systems are parallel one to the other. In particular, when a reference system {t, x, y, z} moves relative 
to another reference system {t,x,y,z} with a velocity parallel to the x-axis, the corresponding 
Lorentz transformation is 



X = /3{x + vt), y — y, z ~ z, t ~ (3{t + vx) , 



(49) 



with the Lorcntz parameter 

Recall that we use a system of units with c — 1. 

Since the different directed Lorentz transformations do not commute, a general transformation 
(with an arbitrary directed vector of velocity) can not be generated by a successive application 
of three orthogonal transformations (relative to three axes). Although, a formula for the Lorentz 
transformation with a general velocity vector is known from the literature, we present here a brief 
derivation from 1)49(1 . This derivation is instructive because the proof conforms with the method 
used in this article. 

Consider a reference system which axes are parallel to the corresponding axes of a rest reference 
system. Let the origin of the reference system move with an arbitrary directed velocity v. Consider 
a radius vector r directed to an arbitrary point in space. Its projection on the direction of v is 

< V r > 

Pvr ^ V ^ '„ ^ . (51) 

V 

Exhibit the vector r as a sum of its tangential and normal parts 

r = Pvr + N^r , iVvr = r - v ^ ^'^ ^ . (52) 

V 

Due to (|49|) . for a Lorentz transformation directed by v, these two parts transform as 

P^r = (3{P^r + vt), (53) 
iVvf = iVvr. (54) 

Consequently, the transform of the spatial coordinates is 

f = Pvf + A^vf = PiP^^r + vi) + N^r , (55) 

or, explicitly, 

f = r + V {pt — a < r,v >) , (56) 

where 

" = 2 — = ~ 1 ^ /, , • (57) 

The change of the time coordinate is also governed only by the tangential part of the vector r 

i ^(3{t+\ \P^r\ \v) =p{t+< P^r, V >) , (58) 

or, explicitly, 

i=P{t+<r,^^>). (59) 

In the special case of a motion parallel to the axis x, the relations (|56|l . (|59|l reduce to the ordinary 
form of the Lorentz transformation H49() . Therefore an arbitrary Lorentz transformation takes the 
form 

*>^(;+</'^>)' , (60) 

r = r + V (/3i - a < r, V >) . ^ ' 

B The linear equation: The 1-point solution moved 

Calculate the d'Alembertian of the field of a single particle moving on a trajectory ■0(i) 

771 

^ = R=\\R\\, (61) 

where 

R = {r - ro) - aij) < ip, (r - rg) > ■ (62) 



p = ^^—^i + lm' + l\m' + --- (63) 



Here a and (3 are Lorentz functions of the velocity n-'-n^ 

1 

/l-ll^ll 

and 

a^-±-ll-^^=]^-Ul-hn'-lUr + --- ] . (64) 



'1 - 

The time derivatives of these functions are correspondingly 



dt 

and 



^-2/3'<tA,i/;>=<t/',t/;> (i + ^\m' + -- - ) , (65) 



^ = 2a'<^,t^>--<^,^> Q + + , (66) 

where a', (3' denote the derivatives relative to the variable 'tp. 
The first order time derivative of the vector H62|l is 

it = - |^2q!'i/) < -i/;,!/? > < '0,r - ro > 

-aip < r - ro > -ip'ip < tp,'4' > -(itp . (67) 
Due to (|63I66|I . up to the order 0{ipip), the second order time derivative of the field is 

§ = '^(<i>,R> + <i>,Rt>)~3^<i>..R><R,R> . (68) 

Thus we get 

= '^<i,.R> + '^{3<i,,R>' -R'\mf)- (69) 

Let us calculate now the spatial derivatives of the field (|61|l . Introduce a set of unit vectors 
61,62,63 directed along the axes x,y,z correspondingly. Thus the x-component of the vector H62I) 
can be written as 

= 61 - at/» < i/), 61 > . (70) 
Consequently, the second order derivative is 

g = -^l|i^.lP + 3^<i^.,i^>^. (71) 

Due to the spherical symmetry of the field ip, its Laplacian takes the form 

= 3-^ (< R^,R>^ + < Ry,R>^ + < R.,R>^) - ^ {\\R,\f + \\Ry\f + \\R.\\^) . (72) 

Since 

< Rx,R > ^ < R, ei > -a < ip , ei >< ip,R> (73) 

we get 

< Rx,R>^ + < Ry,R>^ + < R,.R>^^ -2a< ip,R>^ <ip,R>^ , (74) 



and 



Vy 

Thus the Laplacian of the field ip takes the form 



R,\f + \\Ry\f + \\R,\f = 3 - 2aMf+a^M\^ . (75) 



TYl • • • 

= j^ia'Ml^ - 2a){3 <^Ij,R>^ -WM^R^) ■ (76) 
Correspondingly, the d'Alenibertian of the field is 

□ ^= ^<V',H> +-^(/32 + 2a - o?U\?){i <^,R>' ~R^m') • (77) 
Using the expressions ()63|l for the functions a, /?, we derive 

+ 2a-a^\\ip\\^ ^0. (78) 
Consequently, the d'Alenibertian of the 1-singular ansatz H61|l takes the form 

□^=^<tA,i?>. (79) 

It is clear that a singularity that moves on a straight trajectory with a constant velocity i/? = is 
described by a solution of the linear equation Oip — 0. 

C The non-linear equation: The motion of a 1-point singu- 
larity 

In this section we substitute the ansatz with one singular point — 

1 / 771 \ • • 

ip = --\n[l - k—) , R= {r - r) - atp < ip,{r - r) > -l34> . (80) 



into the left hand side of the non- linear field equation (|33() . From the calculations above, the time 
derivatives of the vector R are 

Rt = -(31P , Rtt = -P-^ ■ (81) 
Thus the second order time derivative is 

As for the spatial derivatives we use 

Rx = ei - aip < ij),ei > , R^^ = , (83) 
where ei is a unit vector along the x axis. Consequently, 

m < R,Rx>^ + < R, Ry>'^ + < R, R^ >2 



Aip = 3 



i?5 1 - fcf 



m_ < R^,R^ > + < Ry, Ry > + < R,,R, > 

km^ < R,R^>^ + < R, Ry>^ + < R, R^ >^ 
(l-fcf )2 

Substituting here the value of Rx we get 

m i?2 - 2a < 1/1,12 >2 +a2||'0||2 < >2 m 3 - 2a 

A(/5 = 3 

k7n^ i?^ - 2a < rp,R>'^ +a^\\xp\\'^ <-ip,R> 
~R6 (1- fcf )2 



(84) 



i?5 i-'tf R^ i-'^f 

2 Z?2 _ / r? \2 1 „,2|U/,||2 / J, jy ^2 



(85) 



Thus 



m(3 S^3<R,^p>'^ +R'^ <R,ip>-PR'^\\ip\\'^ m"^ (3"^ k < R, 
^"P = -f^ T^H + 



i?5 l-'tf (l-fcf)2 

m R^ -2a<ip,R>^ +a'^\\ip\\'^ <ip,R>^ m 3 - 2a\\ip\\^ + a^\\ip\\^ 
R^ l-fcf ^R^ l-yfcf 

/cTO^ i?2 _ 2a < ■ip,R>^ +a^\\i>\\^ <ip,R>^ 

= l^T3fc^ + ]^(^ IHI ) + 

m^fc -i?2 + (/?2 + 2Q-a2||^/>||2)<i^,-0>2 
Using the relation 



(86) 



/3^ = a"||V|r-2a (87) 

we obtain 

771/3 < i?, l/) > 77l2fc 1 

As for the quadratic part of the field equation, to the same accuracy, 



_ 7772/32 < fl^ ^ >2 rft2 i?2 _ 2a < >2 +a2||i/)||2 < I/,,!? >2 

i?6 (i-fcf)2 ~;r^ (1-fcf )2 

7772 1 m? , 0, ,•,,,,< >2 , ^ 

^ -i?J(r3If)I + ^(/5^ + 2a-a^|Hn^^e4)T. (89) 

Using once more the relation H87|) we get 

ab ™^ 1 

When 188(1 and (|90|l are substituted into the left hand side of the field equation 133(1 we obtain 

1-^ , „h < R.'ib > ,„ , 

□ y-fcr;°V,ay,fa = ^ ^Jj^ ■ (91) 

Consequently, the field equation 1(33(1 is satisfied only for a singularity that moves on a straight 
trajectory with a constant velocity. 

D The non-linear equation. The A^^-point solution moved 

Take the field of N singular points as a superposition 
1 ^ / 

= -- ^ In (1 - fc-^ j , R, = {v-v{)- < ^„ (r - r,) > -P,^p.^ . (92) 
i=i ^ ' ^ 

Substitute it in the left hand side of the non- linear field equation ((33(1 . Since this ansatz a superpo- 
sition of N independent solutions, the linear part of the field equation is a sum of the expressions 
given in (P?H|l 

_ m,P, < R„ V?,; > 7772fc 1 



Calculate now the non-linear part 
Thus 

TY mi/3i 

i^tf = H T^l™! '^J' ^ ■ ^^^^ 

As for the spatial derivatives 

< Riy,Ri >< Rjy,Rj > + < RizjRi >< Rjz,Rj > ^ • (96) 

Applying the relation 

Rix^ei~aiipi<ipi,ei>, (97) 
the expression in the brackets of (|96|l takes the form 

(•••) = < R,,Rj > -aj(^<ipj,Ri><ipj,Rj >'j -a,(^<->p^,Rj ><ip^,Ri>'j + 

aia-i < ipi, R, >< ipi^tpj >< -ipj^Rj > ) . (98) 

Consequently, the quadratic part of the field equation takes the form 

"dS" "b^ r 

< Ri , Rj > + 



<4>i,Rj ><ihi,Ri > {p'lPj + a, + - Qfiaj <'>pi,ipj > ) (99) 



Observe that, for i ^ j, this expression coincides with (|9U|I . For i ^ j, the expressions in the second 
line of (|99|l are not canceled, they can, however, be neglected in the lowest approximation. 
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